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ABSTRACT 

We present results of Hubble Space Telescope and Chandra X-ray Observatory observations of globular 
clusters (GCs) and low-mass X-ray binaries (LMXBs) in the central regions of Centaurus A. Out of 440 GC 
candidates we find that 41 host X-ray point sources that are most likely LMXBs. We fit King models to our 
GC candidates in order to measure their structural parameters. We find that GCs that host LMXBs are denser 
and more compact, and have higher encounter rates and concentrations than the GC population as a whole. We 
show that the higher concentrations and masses are a consequence of the dependence of LMXB incidence on 
central density and size plus the general trend for denser GCs to have higher masses and concentrations. We 
conclude that neither concentration nor mass are fundamental variables in determining the presence of LMXBs 
in GCs, and that the more fundamental parameters relate to central density and size. 

Subject headings: galaxies: elliptical and lenticular, cD — globular clusters: general — X-rays: binaries 



1. INTRODUCTION 

Since the discovery of X-ray sources associated with glob- 
ular clusters (GCs) in the Milky Way (MW; e.g., Giacconi 
et al. 1974) it has been known that low-mass X-ray binaries 
(LMXBs) are formed more efficiently in GCs by a factor of 
100 compared to the field (Katz 1975; Clark 1975). This is 
thought to be a direct consequence of the high central den- 
sities in GCs. At high density, dynamical formation mecha- 
nisms are greatly enhanced relative to the field (Clark 1975; 
Fabian et al. 1975; Hills 1976). 

While the dynamical origin of most GC LMXBs was pro- 
posed three decades ago, subsequent observational studies 
of the core properties of GCs that host LMXBs have been 
largely restricted to the MW. The Chandra X-ray Observa- 
tory has made it possible to study the LMXB populations of 
nearby (< 30 Mpc) galaxies (e.g., Fabbiano 2006 and refer- 
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ences therein), but the core properties of the GCs, which are 
thought to be the most relevant in the dynamical processes 
that create LMXB progenitors, are not robustly determined 
even with HST for most of those galaxies. 

In the MW, Bellazzini et al. (1995) showed that bright 
LMXBs appear in denser and more metal-rich GCs, while 
Pooley et al. (2003) and Heinke et al. (2003) showed that the 
number of GC X-ray sources correlates with V, a basic indi- 
cator of dynamical encounter rates defined as V = p\j r^. (e.g., 
Verbunt et al. 2007), where po is the central mass density and 
r c is the core radius. Bregman et al. (2006) showed addition- 
ally that bright GC LMXBs in the MW appear preferentially 
in GCs with smaller core and half-light radii and shorter half- 
mass relaxation times. Sivakoff et al. (2007; S07) used 1 1 
early-type galaxies to show that encounter rates are a good 
predictor for the appearance of LMXBs in GCs, deriving T 
from measurements at the half-light radii. This is consistent 
with earlier results in M87 that used more uncertain derived 
core properties (Jordan et al. 2004; J04). 

To make further progress in understanding the dynamical 
processes that create LMXB progenitors in the cores of GCs 
it is necessary to directly observe LMXBs and the core struc- 
ture of the GCs that host them for as large a sample as possi- 
ble. In this work, we use HST/ACS and Chandra observations 
of the central regions of Centaurus A (or NGC 5128; here- 
after Cen A) in order to extend such studies beyond the Local 
Group. At the distance of Cen A, D « 3.7 Mpc (average of 5 
distance indicators, see §6 in Ferrarese et al 2007), it is possi- 
ble to study the central properties of GCs with HST, an ability 
we exploit in what follows in order to study which structural 
properties of GCs are most relevant in determining the pres- 
ence of LMXBs in GCs. 

2. OBSERVATIONS 

Optical Catalog. We have made use of 21 fields in Cen A 
observed with the F606W filter using the Wide Field Channel 
mode of the Advanced Camera for Surveys (ACS) on board 
the Hubble Space telescope (HST). Nine of these fields were 
observed as part of program GO- 10597 (PI: A. Jordan), which 
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was designed to be combined with previous observations from 
program GO-10260 (PI: W.E. Harris) in order to cover a large 
area in the inner r < 8 kpc plus some fields at larger radii. 

Each field of GO-10597 consists of four 525-sec exposures 
plus a single 58-sec exposure, while the data of GO-10260 
consist of three exposures of 790 sec. Full observational de- 
tails of GO-10260 are given in Harris et al. (2006), but we 
note that for this study we have re-reduced the GO-10260 data 
in parallel with our reductions of the new GO-10597 data to 
construct an independent and homogeneous catalog of GCs. 
A detailed account of the data reduction procedures and GC 
catalog construction will be given elsewhere so only a brief 
summary follows here. 

All data were drizzled using the Apsis package (Blakeslee 
et al. 2003) onto frames with a pixel scale of 0"05. Ob- 
ject detection and photometry were performed using SEx- 
tractor (Bertin & Arnouts 1996). Our detection threshold 
corresponds to »Jf606w ~ 22 AB mag and is such that we 
are not subject to incompleteness effects. We used photo- 
metric zeropoints and extinction coefficients from Sirianni et 
al. (2005) and de-reddened the photometry using E(B-V) = 
0.115 (Schlegel et al. 1998). We dealt with foreground-star 
contamination by eliminating all objects that were consistent 
with the point-spread function (PSF). We matched our object 
catalog with the full photometric catalog (56,674 objects) of 
Peng et al. (2004), obtaining VI photometry for all objects 
for which a match was found. The same catalog was used 
to place our objects on a consistent astrometric frame. For 
objects which did not have a match in Peng et al., we fur- 
ther matched to deep VI photometry of the central regions of 
Cen A obtained with the VET (Minniti et al. 2004). 

An azimuthally averaged surface-brightness profile was ob- 
tained for each object using the ELLIPSE task in IRAF. We 
then fit PSF-convolved King (1966) models to each object fol- 
lowing the procedures applied by McLaughlin et al. (2007) to 
the GO-10260 data (see also Barmby et al. 2007). The best-fit 
models were used to infer various global and core parameters 
for every GC. 14 

To infer mass-based cluster quantities we used F6Q6W 
mass-to-light ratios as a function of (V-/) color derived from 
the models of Bruzual & Chariot (2003) assuming an age of 
13 Gyr for all objects 15 . We culled our catalog by reject- 
ing objects with (V -I) < 0.6, (V-/) > 1.55 or concentra- 
tion c > 2.5. We further eliminated objects whose surface- 
brightness profiles spanned less than 0.5 mag arcsec" 2 or that 
were visually deemed to be background objects. The final cat- 
alog consists of 440 objects, 407 of which have VI photometry 
available. We estimate that ^3% of these objects may be con- 
taminating background galaxies (based on an analysis of 21 
control fields observed in the same filter to similar depths, and 
analyzed in the same fashion as our Cen A program fields). 

X-ray Catalog. Prior to 2007, Chandra observed Cen A 
with the ACIS detectors four times (Observations 316, 962, 
2978, and 3965). In 2007, Chandra performed six deep ob- 
servations (~ lOOks each) of Cen A as part of the Cen A 

14 At the distance of Cen A, a typical GC core radius of r c ~ 1 pc cor- 
responds to 0"056, comparable to the half-width at half-maximum of the 
PSF. Holland, Cote, & Hesser (1999) and Harris et al. (2002) have shown 
that estimates of r c from PSF-convolved King-model fitting in this situation 
are unbiased in general, and Barmby et al. (2007) have shown that the fits 
of McLaughlin et al. (2007) to the GO-10260 GCs in Cen A, which include 
~ 25% of the GCs in our sample, give rise to core-parameter correlations that 
are consistent with those found for MW GCs. 

15 For objects without measured (V —I) we used (M/L)p6ogw = 2-2 in solar 
units. 



Very Large Project (CenA-VLP; PI: R. Kraft). In this anal- 
ysis, we report initial results from an X-ray point source list 
constructed using the first four CenA-VLP observations (Ob- 
servations 7797, 7798, 7799, and 7800) plus the previous ob- 
servations listed above. 

The reduction of the Chandra data prior to 2007 is de- 
scribed in Woodley et al. (2007). The reduction of the CenA- 
VLP observations (a total of 373,353 s) will be presented in 
detail in forthcoming papers; here, we only consider a few as- 
pects relevant to point source detection (see Hardcastle et al. 
2007 for an image of the X-ray data). Sources were detected 
in the 0.5-7 keV X-ray image from each observation using 
WAVDETECT with wavelet scales spaced by a factor of \/2 and 
ranging from 1 to 32 pixels, with a source detection threshold 
of 10~ 7 . Sources that were not associated with the X-ray jet 
or lobes in Cen A were used to register the relative astrom- 
etry between the observations. The majority of the detected 
sources are X-ray binaries in Cen A, and in particular the ones 
associated with GCs are expected to be LMXBs. In what fol- 
lows we will use the term LMXBs for X-ray sources associ- 
ated with GCs. 

Source Matching. After determining the astrometric off- 
sets required to bring the astrometry of our list of X-ray 
sources into the same system as the GCs, we carried out 
source matching with a matching radius of l."6, obtaining 41 
matches. With our chosen matching radius we expect < 1 
false matches. The rms difference between the X-ray and op- 
tical positions for matched sources was < 0."5 in both right 
ascension and declination. 

3. RESULTS 

Previous studies have found that LMXBs appear prefer- 
entially in more luminous (massive) and redder (metal-rich) 
GCs (e.g. S07; Kundu et al. 2007; see Fabbiano 2006 for a 
recent review and Minniti et al. 2004 and Woodley et al. 2007 
for studies in Cen A). Metal-rich GCs are w 3 times more 
likely to host a LMXB, although the scatter around this value 
for different galaxies is significant. 

In the leftmost panel of Figure Q] we show the histogram of 
AB magnitude in the F606W band for all GC candidates (un- 
filled histogram) and for those that have an associated X-ray 
point source (filled histogram). The next panel shows equiva- 
lent histograms in (V-I) color for sources with VI photome- 
try. The numbers in each of the panels of this figure show the 
p-value of a two-sample Wilcoxon test between the two distri- 
butions shown. The p-value gives the probability that the two 
distributions have the same mean values. It is clear that there 
is a strong preference for LMXBs to be hosted by more lumi- 
nous (massive) GCs. The color distribution of LMXB hosting 
GCs is slightly more weighted toward redder GCs, although 
the difference between the mean of the two distributions is not 
significant for our sample. 

In the five rightmost panels of FigureQ]we show histograms 
of central mass density (poX core encounter rate (T), King 
(1966) core (r c ) and half-light (r/,) radii and concentration 
c = log(r, / r c ), where r, is the model tidal-radius, for all GCs 
(unfilled histograms) and those that host an X-ray point source 
(filled histograms). In agreement with observations in ellip- 
tical galaxies (J04; S07) and the MW (Bregman et al. 2006), 
we find that the GCs hosting LMXBs have smaller half-light 
radii and higher encounter rates. Our new observations allow 
us to show explicitly for the first time in an elliptical galaxy 
that they also have significantly smaller core radii and higher 
central densities, as is the case in the MW (Bellazzini et al. 
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FIG. 1 . — This figure shows histograms of total apparent F606W AB magnitude (mF606w)> de-reddened (V—I) color, logarithm of the central mass density (po, 
in units of Mq pc~ 3 ) and encounter rate (r = pi King (1966) core radii (r c ), half-light (r/,) radii and concentration (c) for the full sample of GCs (unfilled 
histograms) and for the subsample that coincide with an X-ray source (filled histograms). The numbers in each of the panels of this figure show the p-value of a 
two-sample Wilcoxon test between the two distributions shown. The p-value gives the probability that the two distributions have the same mean values. 



1995; Bregman et al. 2006). Finally, the concentration is also 
found to be higher in LMXB-hosting GCs (but see below). 

It is well known that many of the structural properties of 
GCs in the MW are correlated with one another (e.g., see 
Djorgovski & Meylan 1994; McLaughlin 2000). It is there- 
fore important to ask which of the parameters in Figure Q] 
might be fundamental in determining the presence or absence 
of LMXBs in GCs, and which parameters might just appear 
to be important because they correlate with more physically 
relevant properties. Of particular interest is disentangling the 
degeneracy between central density (and, thus, the encounter 
rate T = Pq 5 ^) and total mass, which arises from the fact 
that massive GCs are denser on average (McLaughlin 2000; 
Jordan et al. 2005). In principle, other quantities of relevance 
to LMXB production — such as, e.g., the fraction of primor- 
dial binaries — might depend on GC mass but not on either of 
po or r on its own, which leaves the implications of Figure Q] 
somewhat unclear. 

We may empirically probe which GC properties influence 
the presence of LMXBs in GCs after taking into account any 
dependences on total mass, by comparing the properties of 
clusters that host LMXBs with the properties of clusters that 
do not host LMXBs but have the same underlying mass dis- 
tribution. To do this, we first estimate the mean dependence 
of each of Y = {po,T,r c ,ri 1} c} on GC mass M for GCs that 
do not host an LMXB using a robust local smoothing of the 
M-Y scatterplot performed with the Lowess method (Cleve- 
land 1979). We then use the functions Y(M) thus estimated 
to predict the expected average of the variable Y when having 
the same mass distribution as the GCs that host LMXBs by 
computing (Y) = N~ l Y^=\ Y(Mi), where Mi are the masses of 
the N = 4l GCs that host LMXBs. 

We show the result of this exercise in Figure [2] where we 
show various distributions of structural quantities for GCs that 
host LMXBs, indicating with an error bar the 99% confidence 
interval for the median of these distributions and with an ar- 
row the expected median for GCs that do not host an LMXB 
and have the same underlying mass distribution. The result 
of a one-sample Wilcoxon test to estimate the probability that 
the median of the LMXB-hosting sample is consistent with 
that expected for GCs that do not host LMXBs is indicated in 
each panel. 

A very interesting result is given by the rightmost panel in 
Figure |2 after the effects of mass are taken into account, the 



average King concentration shows no statistical difference be- 
tween GCs that host LMXBs and those that do not. In the 
MW there is a trend for more massive GCs to be more con- 
centrated, roue hly following 10 e oc M° 4 (McLaughlin 2000). 
A consistent trend is also present in our sample. We have 
shown that once this trend is taken into account concentration 
per se has no significant effect in determining the presence of 
LMXBs. 

Figure |2] also shows that all of po, T, r c , and r/, have 
mean values that are significantly different for GCs that host 
LMXBs, even when considering clusters with the same mass 
distribution. Thus, processes driven by high core densities, 
and presumably related to high stellar encounter rates, play 
a direct role in enhancing the presence of LMXBs in GCs; 
the dependences on core properties and r/, in Figure Q] are not 
purely coincidental side effects of physics or initial conditions 
that scale fundamentally with cluster mass only. 

Is there a role left for GC mass in determining the presence 
of an LMXB after the dependence on po and r c (through T) 
has been taken into account? We can answer this question by 
repeating the exercise above but now comparing the mass dis- 
tribution of GCs that host LMXBs with those of GCs that do 
not host LMXBs and have the same underlying distribution of 
r. We find when doing this that the average mass of GCs that 
host LMXBs is statistically indistinguishable to that of GCs 
that do not (Wilcoxon p-value of 0.98). As there are no re- 
maining differences in the average masses after dependencies 
in structural parameters have been taken into account through 
r, we conclude that the dependence of LMXB incidence on 
mass is a consequence of a more fundamental dependence on 
central density and size. We note also that c shows no signifi- 
cant difference between GCs that host LMXBs and those that 
do not when considering the same underlying T distributions. 

We follow the procedure described in §4 of S07 to deter- 
mine the dependence of the expected number of LMXBs, A,, 
on GC properties using a maximum-likelihood method. In 
particular, we fit for the exponents in the following assumed 
forms for A,: (1) A, ocM /3 10 l5(l/ " /) r^ and (2) A, oc T^IO^-'K 

We find that the values of 6 obtained are not significantly 
different from zero. This is consistent with the fact that the 
(V-/) color distribution of GCs that host LMXBs is not 
significantly different from the global color distribution in 
our sample. When using the functional form (1) we find 
= 1 .4 ± 0.2, S = ± 0.6 and e = -2. 1 ± 0.4, while using form 
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FIG. 2. — This figure shows histograms of the central mass density (po), encounter rate (r = pb s rf), King (1966) core radii (r c ), half-light (r;,) radii and 
concentration (c) for the full the sample of GCs that are associated with an X-ray point source (X-GCs). The error bar in each panel shows a 99% confidence 
interval for the median of each distribution and the arrow indicates an estimate of the expected value of that median for GCs that do not hold LMXBs (nX-GCs) 
and that have the same undelying mass distribution as the GCs shown in the histogram. The number in each panel is the p-value of a one-sample Wilcoxon 
test to test the probability that the sample median is consistent with the value indicated by the arrow. The ratio of the median values of the quantities that show 
significantly different distributions (in the sense X/nX) are: 2.1 (po), 1.93 (r), 0.68 (r c ) and 0.77 (ry,). 



(2) we find 5 = -0.3 ± 0.6 and r) = 0.85 ± 0.12. These re- 
sults are in good agreement with the results of S07 for Virgo 
Cluster ellipticals, where the analysis was restricted to using 
quantities defined at the half-mass radius. 

4. CONCLUSIONS 

We have shown in this Letter that GCs in Cen A that host 
LMXBs have significantly higher central densities, smaller 
sizes, and higher concentrations than the GC population as 
a whole. The higher concentrations are shown to result from 
the dependence of LMXB incidence on central density and 
size, plus the fact that denser (or more massive) GCs are more 
centrally concentrated on average. 

More importantly, we have used our Cen A sample to com- 
pare the properties of LMXB-hosting GCs to those that do not 
but have the same underlying mass or T distribution. We find 
that GCs that have LMXBs have significantly higher central 
densities, encounter rates and smaller sizes than GCs of the 
same mass that do not have LMXBs. Interestingly, we find 
that the mean concentrations are indistinguishable once the 
trend for more massive GCs to have higher concentrations is 
removed. As r/,/V f depends only on c for King models, our 
conclusion agrees with the MW data presented in Bregman 
et al. (2006; their Figure 1). 

We further show that potential GC mass-dependent pro- 
cesses are not fundamental in the formation of GC LMXBs 



(see also S07; Verbunt et al. 2007). It is rather smaller sizes 
and denser cores, or equivalently higher values of the core 
encounter rate T, that are the main drivers. 

Finally, our finding that concentration is not important in 
determining the presence of LMXBs in GCs validates the use 
of half-mass densities and radii to probe the dynamical prop- 
erties of GCs and their connection to the presence of LMXBs. 
This is helpful as core properties are uncertain for GCs in 
galaxies much more distant than Cen A. The detailed depen- 
dence of LMXB incidence on GC core properties will aid in 
further understanding of the interplay of various dynamical 
processes in creating X-ray binaries (e.g., Banerjee & Ghosh 
2006; Ivanova et al. 2007). In particular, we find in agreement 
with previous work (J04; S07) that the dependence of the ex- 
pected number of LMXBs A, on encounter rate, A, oc T 85 , is 
shallower than the naively expected linear behavior. 
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